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ABSTRACT 

We present the frequency resolved energy spectra (FRS) of the low-mass X-ray bi- 
nary dipper XB 1323—619 during persistent emission in four different frequency bands 
using an archival XMM-Newton observation. FRS method helps to probe the in- 
ner zones of an accretion disk. XB 1323—619 is an atoll source and a type-I burster. 
We find that the FRS is well described by a single blackbody component with kT 
in a range 1.0-1.4 keV responsible for the source variability in the frequency ranges 
0.002-0.04 Hz, and 0.07-0.3 Hz. We attribute this component to the accretion disk and 
possibly emission from an existing boundary layer supported by radiation pressure. 
The appearance of the blackbody component in the lower frequency ranges and disap- 
pearance towards the higher frequencies suggests that it may also be a disk-blackbody 
emission. We detect a different form of FRS for the higher frequency ranges 0.9-6 Hz 
and 8-30 Hz which is modeled best with a power-law and a Gaussian emission line 
at 6.41q3 keV with an equivalent width of 1.615 2 keV and 1.3to'g keV for the two 
frequency ranges, respectively. This iron fluorescence line detected in the higher fre- 
quency ranges of spectra shows the existence of reflection in this system within the 
inner disk regions. The conventional spectrum of the source also shows a weak broad 
emission line around 6.6 keV (Boirin et al. 2005). In addition, we find that the 0.9-6 Hz 
frequency band shows two QPO peaks at 1.4+q ° Hz and 2.8^02 Hz at about 2.8-3.1 a 
confidence level. These are consistent with the previously detected ~ 1 Hz QPO from 
this source (Jonker et al. 1999). We believe they relate to the reflection phenomenon. 
The emission from the reflection region, being a variable spectral component in this 
system, originates from the inner regions of the disk with a maximum size of 4.7x 10 9 
cm and a minimum size of 1.6 x 10 8 cm calculated using light travel time considerations 
and our frequency resolved spectra. 

Key words: accretion, accretion disks-methods:data analysis-stars:binaries:general~ 
stars:neutron-stars:individual:XB 1323—619 -X-rays:general 



1 INTRODUCTION 

XB 1323—619 is an X-ray burster that shows varying inten- 
sity dips that repeat with the orbital period. The source is 
discovered by Uhuru and Ariel V (Forman et al. 1978; War- 
wick et al. 1981) and the dips and bursts are recovered with 
EXOSAT (van der Klis et al. 1985; Parmar et al. 1989). The 
typical dips last about 30% of the orbital cycle, and the 0.3- 
12.0 keV intensity varies creating sharp deep features in time 
with about ~70% decrease in the persistent emission count 
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rate. The source is viewed at i — 60-80° (Frank et al. 1987). 
The dip recurrence interval is not precisely known, with the 
best measurement of 2.94(2) hr (Balucinska-Church et al. 
1999). The source is found to have persistent ~1 Hz quasi- 
periodic oscillations (QPOs) (Jonker et al. 1999). 

A BeppoSAX observation of XB 1323—619 between 
1.0-150 keV reveals persistent emission of the source mod- 
eled by a cutoff power-law with a = 1.48 ± 0.01 and E cu t 
- 44.1 1^ keV together with a blackbody with kT = 
1.77 ± 0.25 keV (Balucinska-Church et al. 1999). A more 
recent INTEGRAL time-averaged 4.0-200.0 keV spectrum 
(combined JEM-X and ISGRI spectra) of XB 1323-619 
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can be best fitted using a blackbody emission of kT — 
1.7 ± 0.3 keV and a Comptonized plasma emission model 
(CompTT, Titarchuck 1994) of kT = 196.4 ± 28.5 keV with 
a seed in photon temperature of 0.4 ±0.2 keV, and a Comp- 
ton optical depth to scattering of r — 0.002-0.007 (Balman 
2009). An RXTE observation below 20 keV reveals a cut- 
off power-law with a — 1.23 ± 0.07 and a blackbody with 
kT — 1.36 ± 0.06 keV for the persistent emission (Barnard 
et al. 2001). A recent Suzaku observation between 0.8 and 
70 keV has a non-dip spectrum consistent with a blackbody 
model of kT = 1.35 ± 0.36 keV and a power-law with an 
a = 1.67 ± 0.1 (Balucinska-Church et al. 2009). Boirin et 
al. (2005) detect Fe XXV and Fe XXVI absorption features 
from this source using the same archival observation ana- 
lyzed in this paper. They find changes in the properties of 
the Fe XXV and Fe XXVI absorption features from per- 
sistent to dipping intervals indicating the presence of less- 
ionized material in the line of sight during dips (ionization 
parameter £ decreases from log(£) of 3.9 ±0.1 to log(£) of 
3.13 ±0.07). It is suggested that absorption lines in dipping 
low-mass X-ray binaries are a result of a circularly sym- 
metric photoionized region on the outer disk and as the 
source is viewed from the absorber at all times, the vari- 
ations in the ionization parameter of this region causes the 
different levels of dipping in these systems (see also Diaz- 
Trigo et al. 2006, 2009). Church et al. (2005) has studied 
the same data set as Boirin et al. (2005) using absorption 
line ratios of Fe XXV/Fe XXVI and curve of growth analysis 
leading to the conclusion that a collisionally ionized plasma 
with kT=31 keV can reproduce the absorption lines in the 
XMM-Newton Spectra. As they find that this temperature 
is close to the accretion disk corona (ADC) temperature 
(they have calculated using the BeppoSax spectrum), they 
propose that the absorption lines in dipping low-mass X-ray 
binaries are produced in the ADC and there is no need to 
invoke a separate location where a highly ionized absorber 
exists. 

In this paper, we aim at deriving frequency resolved 
spectra of XB 1323—619 in an attempt to reveal different 
emission regions and absorbing regions in the system ex- 
tracting the variable component of the emission spectrum 
as opposed to the non-variable component by the described 
method in the paper. We expect that this method should, 
also, reveal the location of the emitting and absorbing re- 
gions in similar (dipping) LMXBs. 



2 DATA AND THE OBSERVATION 

The XMM-Newton Observatory (Jansen et al. 2001) has 
three 1500 cm 2 X-ray telescopes each with an European Pho- 
ton Imaging Camera (EPIC) at the focus; two of which have 
MOS CCDs (Turner et al. 2001) and the last one uses pn 
CCDs (Striider et al. 2001) for data recording. Also, there 
are two Reflection Grating Spectrometers (RGS, den Herder 
et al. 2001) located behind two of the telescopes. We use an 
archival XMM-Newton observation of XB 1323—619 with a 
duration of 50 ksec between 2003 January 29 09:05 UTC and 
2003 January 29 28:58 UTC. A thin optical blocking filter 
was used with the EPIC Cameras. EPIC pn and MOS1 were 
operated in the timing mode whereas MOS2 was operated in 
the full window imaging mode. For the FRS calculations in 



our work, we used only the timing mode EPIC pn data which 
provides a time resolution of 30 usee together with the high- 
est sensitivity (i.e., count rate) and energy resolution among 
the EPIC CCDs. The EPIC pn timing mode is a mitigation 
for pile-up in high count rate sources and allow for pile-up 
free data for rates below 1500 c s _1 compressing the data 
to a single one-dimensional row of length 4'. 4. We analyzed 
the pipeline-processed EPIC pn timing mode data using Sci- 
ence Analysis Software (SAS) version 8.0.5. Data for analy- 
sis (single- and double-pixel events, i.e., patterns 0-4 with 
Flag=0 option) were extracted from a rectangular region of 
87" wide column for the source and the background events 
were extracted whenever necessary from a source free zone 
normalized to the source extraction area (ie.e., same rectan- 
gular area). We calculated a conventional source spectrum 
using the SAS task ESPECGET which automatically creates 
a source spectrum, background spectrum, the ancillary re- 
sponse file and the energy response matrix. These files were 
later utilized to calculate the response files for the frequency 
resolved spectra as the spectra were derived with the method 
explained in the next section. In general, background sub- 
traction was applied for the FRS since EPIC pn background 
contribution at frequencies 0.002-300 Hz would be negligi- 
ble. In the analysis, only persistent emission (as opposed to 
dipping parts of the light curve) were used to calculate light 
curves for the FRS analysis. In order to filter the persistent 
emission user-defined good time intervals were created using 
thresholds on count rates resulting in an effective exposure 
of about 40 ksec. Data were also checked for flaring episodes 
(of the background) and type-I bursts of the source and all 
such occurrences were cleaned with the above procedure. We 
obtained an average count rate of 27.3±0.02 for the source 
with a maximum variation of 25% for the persistent emis- 
sion. This was later used to calibrate systematic errors of 
13% on the FRS spectral bins. 



3 ANALYSIS AND RESULTS 
3.1 Frequency Resolved Spectra 

In this work, the frequency resolved energy spectra were cal- 
culated following the prescription of Revnivtsev et al. (1999) 
and Gilfanov et al. (2003) where they applied this method to 
RXTE data of LMXBs. One assumes that the inverse Fourier 
Transform, the light curve, is a function of an energy depen- 
dent term and a time dependent term and the light curves 
at different energies are related by linear transformations. 
One allows for two components in the source emission spec- 
trum, a constant non-variable part of the source emission 
and flux variations of the second variable component. Spec- 
trum of the variable component is the frequency resolved 
spectrum, FRS. In order to calculate the FRS, we created 
power spectra (PSD-Power Spectral Density) in a certain 
range of energy channels (100-4095) of the XMM-Newton 
EPIC pn data with a pre-determined binning of 50, 100 or 
200 energy channels in the given range for 320-16 seconds 
time segments for each PSD. In each set of four frequency 
ranges, spectra were produced with the same channel bin- 
ning and the duration for the PSD segments. Fast Fourier 
Transforms (FFT) were computed to derive the PSD adopt- 
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ing the normalization of Miyamoto et al. (1991) 
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. In this prescription, t m is the time label for each time bin, 
x m is the number of counts in these bins, N p h is the total 
number of photons in each light curve and C is the average 
count rate in each time segment used to construct PSD. The 
Miyamoto normalization is in units of (rms/mean) 2 /Hz mea- 
suring the variability amplitudes in different energies and 
can be integrated over different frequency range of inter- 
est. The top panel of Figure 1 shows the average PSD of 
the persistent emission light curve for the X M M-N ewton 
observation of XB 1323—619 computed using the Miyamoto 
normalization. Finally, for integrated frequency ranges of in- 
terest, the frequency dependent spectra can be constructed 
by 



F{E k ,f m ) = C k ^P k J^)Af j 

where F is the count rate of the spectrum on the frequency 
f m in the energy channel Et- 

We calculated FRS with the above method in four dif- 
ferent frequency bands of 0.002-0.04 Hz, 0.07-0.3 Hz, 0.9-6 
Hz, and 8-30 Hz. The frequency bands were selected accord- 
ing to the simple criteria that there is good statistics in the 
frequency range and the FRS produced have comparable 
count rates. We also used frequency ranges that complies 
with the different components of the PSD. Among these 
bands, 0.9-6 Hz frequency range (the two narrow Lorentzian 
components of the PSD as discussed below in the next sec- 
tion) includes the frequency of 1 Hz-QPO detected earlier 
from this system (Jonker et al. 1999). 

3.2 The Average Power Spectrum and the 
Detected QPOs 

We averaged several power spectra to create a PSD of 
XB 1323—619 using the Miyamoto normalization in units of 
(rms/mean) 2 /Hz. Figure 1 shows a fitted PSD with a com- 
posite model of three Lorentzians and a power-law following 
the prescriptions of Belloni, Psaltis, & van der Klis (2002), 
yielding a xt of 1.0 for 10 degrees of freedom (the white noise 
level is already subtracted). The bottom panel of Figure 1 
displays the same PSD fitted by a single broad Lorentzian 
producing an unsuccessful fit with a Xv value of 2.9 for 17 
degrees of freedom. The very broad Lorentzian models the 
low frequencies. One of the narrow Lorentzians have a peak 
at 1.4lJ"j Hz with a width of about 0.001 Hz. This is simi- 
lar to the detected 1 Hz QPO from this system. The other 
narrow Lorentzian is at 2.8^02 with a width of 0.003 Hz. 
The significance of the detected Lorentzians can be calcu- 
lated following (Van der Klis, 1989 ; see also Boirin et al. 
2000). We assume that SNR (in a) = P m -P re/ /Ep where 
P m is the power at Lorentzian peak of interest, Pre/ is the 
power of the continuum noise at the same peak frequency, 
and Ep is the error of the frequency bin at the Lorentzian 
peak. We note that the standard deviation of the average 
of the powers at the Lorentzian peak is taken as the error 
since it is larger and more representative rather than a the- 
oretical estimation including power at the peak frequency 
and square root of the multiplication of the number of PSD 
averaged together and the number of raw frequency bins av- 
eraged together. We find that the significance of the 1.4l ' 2 
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Figure 1. The average power spectrum (PSD) of XB 1323— 619 
using a Miyamoto normalization in units of (rms/mean) 2 /Hz. A 
certain white noise level is subtracted. The top panel is the PSD 
of the persistent emission data fitted with a composite model 
of three Lorentzians and a power-law model. The bottom panel 
shows the fit to the same PSD with a single Lorentzian displaying 
a high frequency excess. 



Hz Lorentzian is at 2.8 a and the 2.8_ ' 2 Hz Lorentzian is at 
3.1 a confidence level, respectively. In addition, we calculate 
an integrated rms variability of 15% for the QPO at 1.4 Hz 
and 11% for the QPO at 2.8 Hz. Finally, the continuum of 
the high frequencies is best modeled using a power-law with 
an index ot -QA+° 2 2 . 



3.3 Results of the FRS Analysis 

The derived spectra (FRS) were fitted with a set of mod- 
els comprising a combination of, or solely by, a power-law, 
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a blackbody and a Gaussian emission line models, all ab- 
sorbed by neutral matter with column density Nh (these 
are POWER, BBODY, GAUSS, and WABS models within 
XSPEC, respectively). The fit parameters are displayed in Ta- 
ble 1. For the spectral analysis of the derived FRS, XSPEC 
version 12.5.1 has been used (Arnaud 1996). Spectral un- 
certainties are given at 90% confidence level (A% 2 = 2.71 for 
a single parameter). We grouped the FRS spectral energy 
channels in groups of 2-3 to improve the statistical quality 
of the spectra. During the entire fitting procedure a constant 
systematic error of 13% was used ( see sec. 2). The fits were 
conducted within 0.5-9.0 keV range. Figure 2 shows the fit- 
ted spectra in the four frequency bands. The conventional 
spectral analysis can be found in Boirin et al. (2005) and 
Church ct al. (2005). 

The results show that the best fitting model in the 
0.002-0.04 Hz, and 0.07-0.3 Hz ranges, is a single black- 
body (see Figure 2 top panels) with a Xv value of 1.0 
(see Table 1; The best fit results for each low frequency 
FRS are highlighted in bold-face) . The blackbody tem- 
perature is in a range 1.0-1.4 keV (90% confidence level 
errors). For the two low frequency ranges, the power-law 
(POWER) model or a power-law plus a Gaussian emission 
line (POWER+GAUSS) model fits yield X l values in ex- 
cess of 4.0 for very similar degrees of freedom of 10-13. The 
fits with blackbody and a power-law (BBODY+POWER) 
or a combination of blackbody, power-law and a Gaussian 
emission line (BBODY+POWER+GAUSS) models yield 
similar Xv values to blackbody (BBODY) fits. The added 
POWER model gives a power-law index which is unreal- 
istically steep. Also, the included GAUSS model gives a 
normalization for the line of almost zero (see Table 1). 
Thus, we tested the significance of adding the power-law 
(POWER) and Gaussian emission line (GAUSS) models to 
the blackbody (BBODY) model using FTEST. The FTEST 
probability yields a value in a range of 0.40-0.46 for includ- 
ing the POWER and POWER+GAUSS models along with 
the blackbody model which is too high. As a result, we find 
that the BBODY+POWER or BBODY+POWER+GAUSS 
composite-model fits are redundant. 

The other two FRS in the 0.9-6 Hz, and 8-30 Hz ranges 
were found to be consistent with each other and the fits 
were conducted simultaneously to increase the degrees of 
freedom. We find that the two FRS can be fitted best using 
a combination of a power-law and a Gaussian emission line 
(POWER+GAUSS) model with a xl of 1.1 for 8 degrees of 
freedom compared with any other model or combination of 
models (see Table 1; The best fit results for each high fre- 
quency FRS is highlighted in bold-face). An emission line at 
6.1-6.6 keV (90% confidence level; best fit value is 6.4 keV) 
has been found necessary to improve the power-law fits (see 
Figure 2 bottom panels) . The bottom right-hand panel in 
Figure 2 shows the unfolded photon spectra and the line is 
clearly visible. This is the iron K a fluorescence line indicator 
of reflection and/or reprocessing in the system dominating 
the variability at these frequencies. The width (a) of the iron 
K a line is <0.5 keV for the 0.9-6 Hz and 0.7^^ keV for the 
8-30 Hz ranges. The equivalent width of the line is large and 
it is 1.6t?;2 ke V for the 0.9-6 Hz and 1.3±g;^ keV for the 
8-30 Hz frequency bands, respectively. The fits with a single 
blackbody model, a single power-law model or a combina- 
tion of the two (POWER+BBODY) yield unacceptable fits 



(see Table 1). On the other hand, a blackbody, power-law 
and Gaussian emission line (BBODY+POWER+GAUSS) 
composite model fit gives the same xl with the fitted 
POWER+GAUSS model (see Table 1). The blackbody tem- 
perature is 1.4-2.0 for this fit and the included blackbody 
model has a very low normalization in the spectra. We 
tested the significance of adding the BBODY model to the 
POWER+GAUSS model fit using FTEST. The FTEST proba- 
bility yields a value of 0.42 for including the BBODY model 
along with the POWER+GAUSS model in the fitting pro- 
cedure which is too high. Thus, a fit with a blackbody plus 
a power-law and a Gaussian emission line is redundant. 



3.3.1 Statistical Significance of the 6.4 keV Fe Line 

We carried out a detailed test of significance of existence of 
the Gaussian emission line using Monte Carlo simulations 
(e.g., Porquet et al. 2004; Minuitti & Fabian 2006). For our 
null, hypothesis, we assumed that the spectrum is simply an 
absorbed power-law continuum with the same parameters as 
the absorbed power-law model fitted to real data. We used 
the XSPEC FAKEIT command to create 1000 fake EPIC pn 
spectra corresponding to this model with photon statistics 
expected from ~ 40 ksec exposure and grouped the spectra 
exactly as we have grouped the real spectra. Next, we fit- 
ted each fake spectrum using an absorbed power-law model 
with all parameters let to vary to obtain a \ 2 value. We, 
then added a line, restricting the line with a between 0-1.0 
keV in steps of 0.1 keV and, constraining the line energy 
between 6.0-7.5 keV with a step of 0.2 keV between each 
energy value and fitted the fake spectra, this time, with a 
composite model of a power-law plus a Gaussian line model 
(GAUSS+POWER). This resulted in 1000 simulations of 
X 2 values from each hypothesis. Then, the two x 2 values 
from each hypothesis were subtracted. The resulting distri- 
bution of Ax 2 against cumulative number of occurrences 
is displayed in Figure 3. The curve is a Gaussian centered 
on Ax 2 ~ 34.2 and varies between 25 and 45. The improve- 
ment in the x 2 value as calculated using the (POWER and 
POWER+GAUSS) model fits in Table 1 is 33 for 3 degrees 
of freedom. The probability that one gets a detection of a 
(false) Gaussian line in the fake data which gives an im- 
provement in x 2 equal or better than that obtained in the 
real data is measured by integrating (counting the number 
of occurances) Figure 3 for Ax 2 ^ 33. This integration yields 
a probability of 0.002 implying that the line detection is over 
99% confidence level and that the detection in our FRS is 
not by chance. 



4 DISCUSSION 

FRS technique is applied to Galactic black hole candidates 
(Revnivtsev et al. 1999; Gilfanov et al. 2003; Sobolewska 
& Zycki 2006; Reig et al. 2006) to neutron star LMXBs 
(Gilfanov & Revnivtsev 2005, Revnivtsev & Gilfanov 2006; 
Shrader, Reig & Kazanas 2007) and also to active galactic 
nuclei (Papadakis et al. 2005, 2006). The systematic applica- 
tion of this technique to Z sources in several different states 
of the sources yields soft and hard variable thermal compo- 
nents. The normalization of these components change with 
accretion rate and thus state of the source. The softer of 
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Figure 2. Frequency resolved energy spectra (FRS) of XB 1323—619 in four different frequency bands. The top two panels show the 
FRS fitted with a neural hydrogen absorption model (Galactic) and a blackbody model WABSxBBODY (XSPEC models) of emission in the 
0.002-0.04 Hz (the curve in black color), and 0.07-0.3 Hz (the curve in red color) frequency bands. The top left-hand panel is the count 
rate spectra and the top right-hand panel is the unfolded photon spectra in E 2 f (E) . The small panel under the count rate spectra shows 
the residuals in sigmas. The two bottom panels show the FRS fitted with a neural hydrogen absorption model (Galactic) and a power 
law with a Gaussian emission line model WABSx (POWER+GAUSS) (XSPEC models) of emission in the 0.9-6 Hz (the bottom curve in red 
color) and 8-30 Hz (the top curve in black color) frequency bands. The bottom left-hand panel is the count rate spectra and the bottom 
right-hand panel is the unfolded photon spectra in E 2 f (E). The small panel under the count rate spectra shows the residuals in sigmas. 



the two blackbodies in question is identified with the ac- 
cretion disk and the harder one with the boundary layer. 
Particularly QPOs from such systems at the horizontal and 
normal branch indicates a particular hard FRS revealing 
the expectations that the QPOs originate from the bound- 
ary layer. The temperature of the harder disk component 
is found to change slowly with luminosity, but the bound- 
ary layer temperature remains constant with frequency con- 
firming that the boundary layer is dominated with radia- 
tion pressure. These components appear mostly in the rel- 
atively higher Fourier frequency ranges of spectra whereas 
the disk-blackbody components appear in the lower Fourier 
frequency ranges. 

A systematic FRS analysis of an atoll source (4U 1728- 
34) performed by Shrader et al. (2007) shows that the island 
state of the source is dominated by a power law or a black- 
body (kT = 1.4-1.7) plus a power law model that does not 
vary with Fourier frequency for 0.008-0.8 Hz, 0.8-8 Hz, and 
8-64 Hz frequency ranges. They do not detect a resonance 
iron line around 6.4 keV in the FRS. The higher accretion 
rate states of the source at the lower banana and upper 



banana show only a harder blackbody with a temperature 
around kT = 2.0-2.2 keV indicating that at higher luminosi- 
ties and accretion rates the variable spectral component is 
coming from the boundary layer only and is fairly constant 
in normalization and temperature across frequencies and dif- 
ferent luminosities. Gilfanov et al. (2003) study the FRS of 
the atoll source 4U1608-52 in the frequency ranges corre- 
sponding to the frequency average spectrum, kHz QPOs and 
a 45 Hz QPO spectra. Their results show that the FRS is 
consistent with two components one of which is the emission 
from the boundary layer with kT=2.0-2.2 keV or a CompTT 
model of Comptonized plasma emission. However, they find 
that the QPO spectra are mainly consistent with 2.1-2.4 
keV blackbody spectra revealing the radiation pressure dom- 
inated boundary layer emission. 

Our results indicate a blackbody model of emission 
with kT=1.0-1.4 keV for the two lowest frequency ranges 
from 0.002 to 0.3 Hz. This is consistent with the findings 
of Shrader et al. (2007) in similar frequency ranges with 
this study. The conventional spectral analysis of the XMM- 
Newton data shows the existence of a blackbody model of 
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Table 1. Best-fit spectral parameters of FRS for the four different frequency ranges obtained using a power-law, a BBODY model or the 
combination of the two together with a neutral hydrogen absorption model WABS. Also, an additive Gaussian model is used. Nu is the 
absorbing column in units of 10 22 atom cm -2 , kT is the blackbody temperature in keV, K B is the normalization of the blackbody 
model, r is the power law photon index, K pi is the model normalization photon cm -2 s -1 keV -1 at 1 keV, Eq is the central energy 
of the Gaussian model in keV, a is the width of the Gaussian line in keV, Kq is the normalization of the Gaussian line in total 
photons enr 2 s~ 1 keV -1 . The best fitting model to the data in each frequency range is highlighted in bold-face. 
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emission as the softer X-ray component with kT=0.92-1.2 
keV (Boirin et al. 2005; Church et al. 2005) which is, then, 
a variable component in the spectrum of the source. How- 
ever, we note that any FRS calculated in this work shows 
-Ffrs < 0.12F^ where F x is the total X-ray flux inferred 
from conventional spectral analysis. The X-ray luminosity 
of XB 1323-619 (5.2xl0 36 erg s _1 ; Boirin et al. 2005) cor- 
responds to a mass transfer rate of 0.08M e dd (The distance is 
10-20 kpc and M could be as high as 0.13M e dd) • We suggest 
that this soft component originates in the boundary layer 
and is persistent over the two frequency ranges supported 



by radiation pressure as expected. In addition, INTEGRAL 
(an average spectrum over two years, Balman 2009), RXTE 
(Barnard et al. 2001) and Suzaku (Balucinska-Church et al. 
2009) spectral analysis indicate similar blackbody temper- 
atures between 1.0-1.7 keV (overlapping within 95% confi- 
dence level error range) taken at different epochs and rela- 
tively different luminosities in a range 1-5.2 xlO 36 erg s _1 . 
This supports the stability of an existing boundary layer 
emission (see also sec. 1). On the other hand, we need to 
stress that the appearance of blackbody emission in the low 
frequencies such as in this case is a typical feature of the 
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FRS of Z sources (e.g., Gilfanov & Revnivtsev 2005) and 
the blackbody is, then, a disk-blackbody. However, a disk- 
blackbody with such temperatures would be atypically high 
for an Atoll source (like XB 1323— 619) (e.g., Tarana, Baz- 
zano & Ubertini 2008; Di Salvo et al. 2000). 

We find a different form of FRS and hardening of the 
spectra in the two higher frequency ranges. The existence 
of a power-law spectra with an iron fluorescence K a line at 
6.4l ' 3 keV in the two frequency ranges 0.9-6 Hz and 8-30 
Hz with the lowest normalization reveals reflection from a 
media located in the vicinity of the main source of emission. 
In the case for reflection from optically thick cold neutral 
medium, the main reflection features are a narrow unshifted 
iron fluorescence K a line at 6.4 keV, an absorption edge 
at 7.1 keV and a reflected continuum emission peaked at 
20-30 keV (Basko et al. 1974; George & Fabian 1991). How- 
ever, deviations from this reflection model (e.g., broad lines 
and smeared edges) is detected in X-ray binaries indicating 
complicated ionization stages in or motion of the reflecting 
media (e.g., Miller 2007, Cackett et al. 2009a, 2009b, 2008; 
Reis, Fabian & Young 2009; di Salvo et al. 2009; Reynolds 
et al. 2010). The complicated modelling of these modified 
lines show a line centered mainly on 6.4 keV (errors in a 
range from 6.0 to 6.9 keV) with large broadening of the 
widths (equivalent widths in a range 50-400 eV) including 
extended red-wings due to general relativistic effects (i.e., 
redshifts) or Comptonization in a hot inner corona. In the 
case for XB 1323—619 the spectral parameters on Table 1 
reveal that a cold/warm reflection model with the unshifted 
iron Ka line is plausible and in the 0.9-6 Hz range. However, 
we calculate a a < 0.5 keV for the width of the iron fluores- 
cence line in the 0.9-6 Hz range and a a of 0.7to's keV in the 
8-30 Hz frequency band at a 90% confidence level. The error 
on the central energy of the line (i.e., 6.1-6.6 keV) and the 
large width of the line (particularly in the 8-30 Hz range) 
allows for shifted and broad resonance Fe line in the FRS, as 
well. We detect broadening of the line width with increasing 
frequency as was detected in the FRS analysis of Cyg X-l 
by Revnivtsev et al. (1999). In general, the parameters of 
the line is similar to the modified 6.4 keV Fe fluorescence 
lines detected in other X-ray binaries. 

We calculate that the neutral hydrogen absorption, also, 
increases by a factor of 3-4 in the spectra at the two high fre- 
quency ranges. This strongly suggests that the FRS method 
may have revealed the location of an absorbing medium in 
the persistent emission. Particularly, the 0.9-6 Hz FRS shows 
the highest absorption due to neutral Hydrogen as modeled 
in Table 1, but since the FRS have low spectral resolution, 
the true absorption may be of any kind, warm or cold. 

An emission line at 6.6 keV is also detected in the 
conventional spectral analysis of the persistent spectrum 
(Boirin et al. 2005). The XMM-Newton analysis reveals a 
weak broad line with a FWHM of 2.0t°;4 keV - As mentioned 
in the above paragraphs processes like relativistic broad- 
ening, Compton Scattering, rotational velocity broadening 
may be the reason for this large line broadening. The maxi- 
mum limits of the line width we calculate from the FRS are 
smaller than the FWHM derived from conventional spectral 
analysis. Since the line in the conventional spectrum is weak, 
it is sensitive to the choice of the continuum emission. How- 
ever, the large FWHM is a strong indication that it may 
be a blending between the iron fluorescence lines originat- 



ing from the different parts of the disk (i.e., inner and outer 
disk). 

We believe that the two high frequency band spectra 
belong to the reflected emission and thus show a different 
spectra in comparison with the rest of the frequency bands 
that show blackbody emission from the disk. The conven- 
tional spectral analysis of the XMM-Newton data shows a 
power-law emission component as the harder X-ray compo- 
nent with r=l. 6-2.0 (Boirin et al. 2005; Church et al. 2005) 
which is consistent with our results. 

We find that the PSD of the persistent emission shows 
two Lorentzian components around 1.4 and 2.8 Hz at about 
99% confidence level. We belive this is consistent with the 
~ 1 Hz QPO previously detected from this system (Jonker 
et al. 1999). There has been two suggestions for the ori- 
gin of the 1 Hz QPO: 1) quasi-periodic obscuration of the 
central source by a structure in the disk producing rela- 
tively energy- independent oscillations (Jonker et al. 1999); 
2) global disk oscillations (Titarchuck &: Osherovich 2000). 
Jonker et al. (1999) suggest that the QPO can not be ex- 
plained with Compton scattering alone since that would re- 
sult in strong energy dependence of the QPO which is incon- 
sistent with their findings. They suggest that the relatively 
low energy dependence may arise in a model with an in- 
termediate temperature structure containing hot and cold 
electrons combining Compton scattering with absorption. 
Quasi-periodic obscuration of the central source emission by 
an opaque medium, an orbiting bump on the surface of the 
disk, may cause the QPO. The QPO exists in the persistent 
emission, the dips and the bursts consistently revealing the 
connection to our work that the source emission is incident 
on a reflector carrying its signature. The high inclination 
angle of the system also plays an important role. The high 
frequency FRS (0.9-30 Hz) derived in this study supports 
this scenario for the QPOs, since FRS is best fitted with a 
power-law confirming the Comptonizing nature and the ab- 
sorption is then supported by the existence of the 6.4 keV Fe 
K a emission line and higher N^ parameter values from the 
fits. In the light of the spectral results, the QPOs, and the 
6.4 keV Fe K Q emission line, we suggest that high frequency 
FRS relates directly to reflection/reprocessing and scatter- 
ing in the disk and the frequency bands used to derive the 
FRS would yield a finite light crossing time of the reflector 
Tref—^-ref/c, also indicating the location of this region within 
the disk. For f=l-30 Hz, l re / is c/27rf yielding a maximum 
size of 4.7xl0 9 cm for the reflector and a minimum size of 
1.6 xlO 8 cm that would also denote a minimum observable 
inner disk size (the inner disk is at or closer than this limit). 
We suggest that the FRS method applied in this work have 
recovered the reflected component in the system that had 
not been found before by conventional methods. 

XB 1323-619 exhibits iron Fe XXV and Fe XXVI ab- 
sorption lines and the system is viewed relatively close to 
edge-on. In general, this is true for most of the dipping 
LMXBs. The conventional spectral analysis indicates a lack 
of orbital phase dependence of features except during dips 
suggesting that a highly-ionized absorber can be located in 
a thin cylindrical geometry around the compact object. This 
suggested location of the highly ionized absorber in LMXBs 
is mainly in the outer disk (Boirin et al. 2005; Diaz-Trigo et 
al. 2006). In comparison with this scenario, our FRS analysis 
indicates that there is an absorbing region which is associ- 
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ated with the reflection zone in the inner regions of the disk 
at 4.7-0.2xl0 9 cm. However, we note that though we find an 
absorbing medium in the inner regions of the disk, we do not 
rule out another absorbing region (possibly highly ionized) 
which can exist somewhere in the outer disk. A different ap- 
proach using dip-ingress timing for all the dipping LMXBs 
(Church & Balucinska-Church 2004) indicates existence of 
extended accretion disk coronae (ADCs) having radial ex- 
tend typically 50000 km or 5-50% of the accretion disk radius 
(Balucinska-Church et al. 2009). The calculated size of an 
extended ADC in XB 1323-619 is ~2.7xl0 9 cm (Church & 
Balucinska-Church 2004) with a temperature of about 44 
keV. INTEGRAL results show a temperature of about 196 
keV in 4-200 keV range ( Balman 2009) with an estimated 
location of a suggested ADC at r < 1 x 10 9 cm. The assumed 
location of the ADC have some overlap with the reflection 
region size we calculated using our FRS. Then, the 6.4 keV 
line broadening and may be due to Comptonization in an 
ADC. We note that the size of the reflection region is 4.7- 
0.2 xlO 9 cm as mentioned above and the general relativistic 
effects are detected within 40-1 xlO 7 cm (e.g., Papitto et 
al. 2009; Cackett et al. 2009b; Reis et al. 2009) and may 
not explain the large equivalent widths and extended line 
widths. We note that (using the INTEGRAL results) such 
a hot ADC may possibly show band limited noise effects in 
harder energies and at relatively higher frequencies than the 
ones studied in this work. 

Finally, we stress that we do not recover any absorp- 
tion features (i.e., lines) in the four FRS we have calcu- 
lated in this study which may rule out the existence of a 
highly ionized absorber within the inner zones of the ac- 
cretion disk. However, if there were narrow weak absorption 
lines restricted to width a < 0.14 keV and equivalent widths 
< 0.04 keV (Boirin et al. 2005), our crude spectral resolu- 
tion might have missed them out since we detect a strong 
Fe fluorescence line at 6.4 keV. 



5 SUMMARY AND CONCLUSIONS 

We have analyzed the frequency resolved energy spectra 
(FRS) of the low-mass X-ray binary dipper XB 1323—619 
during persistent emission in four different frequency bands 
using an archival XMM-Newton observation. We recover 
two different continuum shapes for the four frequency ranges 
we have studied. We find that the FRS shows a single black- 
body component with kT in a range 1.0-1.4 keV for the 
variability in the frequency ranges 0.002-0.04 Hz, and 0.07- 
0.3 Hz. We attribute this component to the accretion disk 
and possibly emission from an existing boundary layer sup- 
ported by radiation pressure. On the other hand, we stress 
that the emergence of a blackbody component in the lower 
frequency ranges and disappearance towards the higher fre- 
quencies suggests that it may also be a disk-blackbody emis- 
sion given the FRS properties derived for Z-sources. We find 
a different form of FRS for the frequency ranges 0.9-6 Hz 
and 8-30 Hz which is a power-law model together with a 
Gaussian emission line at 6.41q3 keV with an equivalent 
width of 1.6^12 keV and 1.3to' 9 keV for the two frequency 
ranges, respectively. The width a of the line is < 0.5 keV 
and 0.71q g for the 0.9-6 Hz and 8-30 Hz ranges. In addition, 
we recover different neutral hydrogen column density for the 



a) 
n 




Figure 3. A% 2 against cumulative number of chance of occur- 
rences for adding a single Gaussian line in 1000 randomly sim- 
ulated spectra. The null hypothesis spectra is a power-law with 
a neutral hydrogen absorption (Galactic). The Ax 2 values vary 
between 25 and 45 for 8 degrees of freedom. This line detection 
is over 99% confidence level in the frequency resolved spectra. 



FRS derived in our work. The Ah is highest in the 0.9-6.0 Hz 
range with a value of 2.8t og x 10 22 cm -2 and is only slightly 
lower for the 8.0-30.0 Hz range, but lower by a factor of 3- 
4 in the 0.002-0.04 Hz and 0.07-0.3 Hz ranges. So, there is 
more associated absorption towards the inner disk and a dis- 
tinct absorption difference for the two different continuum 
detected for the FRS. We find that the emission structure 
within the inner regions of the disk is such that there is either 
boundary layer emission and an embedded reflection region 
with higher absorption towards the inner regions or there is 
a reflection region and an ADC in the inner disk followed 
by a disk-blackbody towards the outer disk regions. The 
derived iron fluorescence line shows the existence of reflec- 
tion/reprocessing in this system within the inner disk region. 
The size of the reflection region is in a range 4.7-0.16xl0 9 
cm assuming light travel times. We suggest that the FRS 
method applied in this work have recovered the reflected 
component in the system that had not been found before by 
conventional spectral methods. We stress that FRS can be 
used in this manner to find reflection components in other 
LMXBs where conventional spectra can not be adequate. 
Finally, we find that the 0.9-6 Hz frequency band shows two 
QPO peaks at 1.4±*;° Hz and 2.8j£ 2 , Hz at about 2.8-3.1 
a confidence level with integrated rms variability of 15% 
and 11%, respectively. We relate the origin of the QPOs to 
the reflection phenomenon. The width of the 6.4 keV line 
increases with increasing frequency where the location, the 
QPOs are produced, is the outer most regions of the reflec- 
tion zone with the least spread in the line and the highest 
Ah, thus, possibly a cold reflection zone (i.e., region with 
large clumps). 
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